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ABSTRACT 

Context. Globular cluster are believed to boost the rate of compact binary mergers that may launch a certain type of cosmological 
gamma-ray burst (GRB). Therefore globular clusters appear to be potential sites to search for remnants of such GRBs. 
Aims. The very-high-energy (VHE) gamma-ray source HESS J 1747-248 recently discovered in the direction of the Galactic globular 
cluster Terzan 5 is investigated for being a GRB remnant. 

Methods. Signatures created by the ultrarelativistic outflow, the subrelativistic ejecta and the ionizing radiation of a short GRB are 
estimated for an expected age of such a remnant of t > 10 4 years. 

Results. The kinetic energy of a short GRB could roughly be adequate for powering the VHE source in a hadronic scenario. The age 
of the proposed remnant estimated from its extension possibly agrees with the occurrence of such events in the Galaxy. Subrelativistic 
merger ejecta could shock-heat the ambient medium. 

Conclusions. Further VHE observations can probe for the presence of a break towards lower energies expected for particle acceler- 
ation in ultrarelativistic shocks. Deep X-ray observations would have the potential to examine whether there is any thermal plasma 
heated by the subrelativistic ejecta. The identification of a GRB remnant in our own Galaxy may also help to explore the effect of 
such a highly energetic event on the Earth. 
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1. Introduction 

Gamma-ray bursts (GRBs) are enigmatic explosions, detected 
usually at cosmological distances. They appear to come in 
two distinct flavor s according to their duration (see e.g. 
iGehrels et al.L l200"9l for a review). Powerful explosions usually 
leave long-living remnants behind that potentially could be stud- 
ied in the local Universe as relics of the original events. 

Long bursts (duration longer than about 2 s) are gener- 
ally considered to be launched by the death of a massive star 
and several structures in our Galaxy have been suggested to 
be the remnants of such events. These proposed remnants are 
in many cases linked to sources of very-high-energy (VHE, 
>100 GeV) gamm a-ray emission. The VH E gamma-ray source 
HESS J1303-631 dAharonian et al.L I2005I) has been suggested 
as the remnant of a long GRB that happened > 10 4 years ago 
dAtovan et all I2006I) . Also the str ucture W49B has been sus- 
' pected to be a remnant of a GRB dloka et al.L [2004). Recently 
it has been arg ued that the populatio n of unidentified TeV 



sources (see e.g. lAharonian et all 2008I) may be do minated by 
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GRB/Hypernova remnants ( Ioka & MeszarosLl2010l) . 

For short GRBs (duration shorter than 2 s), compact binary 
mergers as ce ntral engine have been identified as t he preferred 
scenario (e.g. IGehrels et all I2005L iLee et all I2005I) . Such com- 
pact binaries may be efficiently created in globular clusters since 
these environments feature high densities of very old stars in 
their cores. Indeed in the globular cluster M 15 already a close 
syste m that consists of two neutron stars (NS) has been discov- 
ered dAnderson et all 119901) . Based on these facts, it has been 
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argued that a considerable fraction of all NSNS binaries are 
formed in globular clusters dGrindlav et all l2006h . It has been 
claimed that the rate of short bursts in the local Universe is dom- 
inated by merger s of dynamically formed compact binaries i n 
globular clusters dSalvatera et al.L I2008L iGuetta & Stelli 120091) . 
Therefore, globular clusters might be the prime environment for 
seeking remnants of short GRBs. Potential signatures of rem- 
nan ts of compact binary mergers hav e already been discussed 
(see iDomainko & Ruffertl 120051 120081) . 

The H.E.S.S. collaboration has very recently reported the 
detection of a VHE gamma-ray source HESS J1747-248 
in the direction o f the galactic globular cluster Terzan 5 
(lAbramowski et all 1201 lh with a flux above 440 GeV of 
(1.2+0.3)xl0~ 12 cm _2 s _1 . The source is extended with intrin- 
sic extension of 9'. 6+2'. 4, and there appears to be an offset 
from the cluster core position by ab out 4'. 0+1'. 9 . Terza n 5 
is located at a distance of 5.9 kpc dFerraro et all 120091) . at 
RA(J2000) 17 h 48 m 04 s .85andDec -24°46'44".6 (1 3.8°, b 1.7°), 
and exhibits a core radius of 0.15' a ha l f-mas s radius of 0.52' 
and a tidal radius of 4.6' dLanzoni et al.L |2QToT) . Terzan 5 is the 
globular clu ster with the h i ghest expected rates of close stellar 
encounters dPoolev & Hutt 120061) and with the larges t popula- 
tion o f millisecond pulsars (33) discovered up to now dRansomL 
120081) . Gamma rays of likely magnetospheric origin produced 
by millise cond pulsar s have been detected in th e GeV range in 
Terzan 5 dKong et all [20ToL lAbdo et all l2010al) . Diffuse X-ray 
emission extending beyond the half-mass radius has also been 
reported dEger et al.L 120101) . In the vicinity of T erzan 5 several 
struct ures in the radio band have been found dClapson et all 
Terzan 5 is a predicted VHE gamma-ray emitter where 
energetic electrons produced by the large population of millisec- 
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ond pulsars up-scatter stellar photons to gamm a-ray energies 
(IVenter et all l2009t iBednarek & SitarekL 120071) . In contrast to 
these models here VHE gamma-ray emission originating from 
collisions of hadronic cosmic rays with ambient target nuclei 
and subsequent n° decay is explored. A short GRB is adopted 
as the accelerator of the cosmic rays. 

The paper is organized as follows: In Sect.[2]potential indica- 
tions for a hadronic VHE gamma-ray production in Terzan 5 are 
discussed, in Sect.[3]signatures left behind by the ultrarelativistic 
outflow of an ancient GRB in Terzan 5 are assessed, in Sect. [4] 
signatures related to the subrelativistic ejecta are explored, and 
in Sect. [5] potential traces of ionizing radiation emitted by the 
GRB are investigated. 



2. Potential indications for a hadronic scenario 

Some of the properties of the VHE gamma-ray source that were 
detected in the direction of Terzan 5 appear to challenge lep- 
tonic models for the gamma-ray production. In particular, the 
extension, indication for an offset of the source from the glob- 
ular cluster core and a power-law spectrum are not self-evident 
for such a scenario. The intensity of leptonic inverse-Compton 
(IC) radiation scales linearly with the energy density of the target 
photon field. For Terzan 5 the energy density of the stellar pho- 
ton field drops from about 1000 e V/cirT 3 in the core r egion to 40 
eV/cirT 3 at the half-mass radius (IVenter et all l2009h to about a 
few eV/crrT 3 at the extension of the VHE source. Consequently 
a very centrally peaked source centered on the GC would be ex- 
pected, which does not seem to be supported by the H.E.S.S. 
observations. IC emission in the VHE range should be accom- 
panied by synchrotron emission in the X-ray band. Diffuse X- 
ray emission centered on Terzan 5 of possible syn chrotron ori- 
gin has indeed been discovered (lEger et all 1201 Ol) but the po- 
tential offset of the VHE gamma-ray emission from the peaks 
of the X-ray emission and radiation field challenges a leptonic 
scenario. Furthermore, since the optical to near-infrared starlight 
photon field should be up-scattered by the very-high energy elec- 
trons, Klein-Nishina (KN) suppression of the IC process should 
be significant at multiple TeV energies causing a steepening in 
the VHE gamma-ray spectrum. F or a target stellar ph oton field 
with mean temperature of 4500 K (I Venter et all l2009h and elec- 
tron energies of 10 TeV the KN suppre ssion factor is already 
about 0.025 dCoppi & Blandfordl 1 19901) . so the VHE gamma- 
ray spectrum should steepen well before this energy. In compar- 
ison, the observed spectrum may follow a straight power law, but 
this result is influenced by limited statistics. To account for the 
aforementioned arguments, as an alternative to an IC scenario, 
hadronic gamma-ray production is explored in this paper as the 
origin of the VHE source. 



3. Ultrarelativistic outflow 

3.1. Energetics 

GRBs are generally believed to be caused by a pair of ultrarela- 
tivistic jets that are ejected from the central engine. Relativistic 
shock waves accelerat e all particles from the inco ming plasma 
to relativistic energies dBlandford & McKeeL fl976h . thus a sub- 
stantial fraction of the initial energy of the relativi s tic bla st wave 
is transferred into cosmic rays (see lAtovan et all l2006i for the 
case of a GRB remnant). Therefore, for such a scenario, the en- 
ergy in cosmic rays is a measure for the kinetic energy of the 
relativistic outflow. 



From the luminosity of the VHE gamma-ray source the en- 
ergy in hadronic cosmic rays can be estimated if the density of 
target material is assumed. At the location of Terzan 5 the den- 
sity of target materia l should be on the order of n « 0. 1 crrT 3 
dClapson et alll20Tlh . Additionally, to constrain the total energy 
in hadronic cosmic rays for the entire relevant cosmic ray en- 
ergy range above about 1 GeV a spectral index for the region 
below the range that can be probed with H.E.S.S. has to be 
assumed. If a cosmic ray spectral index of 2.0 is adopted be- 
low 5 TeV, the cosmic ray energy that produces gamma rays 
at the energy threshold of 440 GeV of the H.E.S.S. measure- 
ments, the total energy i n hadronic cosmic rays wo uld be £cr ~ 
10 51 («/0.1cm - 3 )- 1 ergs (lAbramowski et al.L[20Tlh . 

Particles accelerated by a relativistic shock wave will 
not feature a single power-law spectrum but should have a 
break towards low energies where the break energy E\x is 
given by the bulk Lorenz factor T of the relativistic shock 
dBlandford & McKeeL Il976t iKatzt 11994 according to E bI ~ 
m p c 2 r 2 /2 with m v the mass of the particle and c the speed of 
light. If a proton energy of 5 TeV, which produces gamma rays 
at the H.E.S.S. threshold, is adopted as the break energy, then 
the Lorenz factor at the time when most particles are accelerated 
would be < 100. Lower break energies would result in lower 
Lorenz factors. A spectral break towards lower particle energies 
would reduce the £cr needed to explain the V HE source de- 
pendi ng on the break energy by up to a factor of 2 dAtovan et all 
120061) . Since ultrarelativistic blast waves are expected to trans- 
fer a substantial part of their kinetic energy into cosmic rays the 
energetics of short GRBs could be roughly adequate for the ob- 
served VHE source if a ratio of prom pt electromag netic energy 
release to kinetic energy of 0.1 - 0.01 dNakarl [20071) is assumed. 
Extending VHE observations t o low er energies with the planned 
CTA array dCTA Consortiuml 1201 Ol) could probe for a break in 
the gamma-ray spectrum of this source to test this scenario. 

3.2. Age of the remnant and rate of short GRBs 

After the acceleration, cosmic rays will diffuse away from 
the location of the GRB and will therefore form extended, 
center-filled gamma-ray sources dAtovan et all 120061) . The age 
of a GRB remnant would in such a case be given by the 
diffusive propagation time of cosmic rays to the extension 
of the source. For the VHE source in Terzan 5 an age of 
the remnant of 10 3 (£)/10 28 cirrs -1 ) -1 years would be found 
dAbramowski et all 1201 ll) . where D is the uncertain diffusion 
coefficient here compared to the value e stimated for 5 TeV p ro- 
tons in the galactic disk of 10 28 crrrs -1 dAtovan et alll2006l) . 

The age obtained for the potential GRB remnant at Terzan 5 
can be compared to the rate of compact binary mergers in the 
Galaxy and to the rate of short GRBs in the local Universe. 
From field NSNS binaries a g alactic merger rate to o ne event per 
(0.5 - 7)xl0 4 years is found dKalogera et all 120041) . For merger 
induced bursts which are formed in globular cluste r s a lo cal 
rate of 20 - 90 events per Gpc~ 3 yr~' (ISalvatera etaD, 120081) or 
~4 Gpc 3 yr 1 dGuetta & Stella! 120091) has been estimated. With 
a density of Milky way- type galaxies in t he local Universe of 
0.01 galaxies per Mpc 3 dCole et all 1200 ll) this results in a rate 
of short bursts per galaxy of about one event per (0. 1 - 0.5)xl0 4 
(/tT'/lOOr 1 years or 2.5xl0 4 (f^/lQQT 1 years. Here / b is the 
beaming fa ctor of short bursts, uncertain in the range of 1 <k 
f b l < 100 dNakarl 120071) . It appears that the age of the poten- 
tial GRB remnant would be roughly comparable to the rate of 
short bursts in the Milky Way if the diffusion coefficient D were 
below 10 27 cm 2 s _1 in the TeV range. Slowed down diffusion at 
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the location of Terzan 5 might be a reasonable assumption since 
X-ray obs ervations demonstrated that the plasma there is hi ghly 
turbulent ( lYao & WangL 120071: ICrocker & Aharoniani 1201 lb . To 
conclude, for possible diffusion coefficients the age estimated 
for the hypothetical remnant could be roughly comparable to the 
merger occurrence. 

3.3. Potential multiwavelength signatures 

Nonthermal diffuse X-ray emission from Terzan 5 ex tending be- 
yond the half-mass radius has also been detected dEger et al.L 
l2010l) . The surface brightness of this X-ray emission roughly fol- 
lows the surface brightness of the stellar component of the glob- 
ular cluster, which could indicate an IC origin. In the framework 
of a GRB remnant scenario these diffuse X-rays might be emit- 
ted by primary electrons accelerated by the ultrarelativistic blast 
wave. Primary electrons would have cooled down below GeV 
energies since their production and th erefore would not emit 
synchrotron radiation in the GHz band (I Atovan et all l2QQ6h but 
might still be energetic enough to scatter the intense stellar radi- 
ation field in the globular cluster to the X-ray range. For such a 
scenario an energy i n electrons of 5 x 10 49 (i/ ra d/40 eV cm 4 ) ergs 
would be required dEger et aTl 1201 Ol) with u m d being the mean 
energy density in the stellar radiation field in the diffuse X- 
ray source. If this situation is indeed realized in Terzan 5, then 
a comparison to the energy in hadronic cosmic rays required 
to explain the VHE gamma-ray emission can be made. Here a 
quite large electron-to-proton ratio of « 0.1 even after the elec- 
trons have cooled down below GeV energies would be found. 
Observations in the hard X-ray regime to extend the spectrum of 
the X-ray source might shed more light on its origin. 



about T * 10 7 K for a mean density of n m 0.1 cm 3 , which 
should radiate 4 x 10 33 (T/10 7 K)~ - 6 («/0.1cm~ 3 ) 2 ergs s" 1 in 
X-rays, corresponding to 10~ 12 (r/10 7 K) ~ 0,6 (n/O.lcirT 3 ) 2 ergs 
cm -2 s _1 at Earth (X-ray emissivity from lSarazi nl ll986h . Such 
an extended and faint structure might have been missed by the 
present X-ray observations due to the limited field of view of 
Chandra and the limited sensitivity of ROSAT. Deep future X- 
ray obse rvations with the ex isting XMM and the planned IXO 
mission dBarcons et all [20TTh would have the potential to probe 
the presence of diffuse thermal X-ray emission in the direction 
of Terzan 5. 

4.2. Signatures from the ejecta 

The ejecta o f compact binary mergers ma y consist of r-process 
nuclei (e.g. lLattimer & Schramm! Il974t iRuffert et ail 11997b 
iFfeiburghaus et all 11999b ." Many of these heavy nuclei are ra- 
dioactive and emit hard X-ray and soft gamma-ray line emission 
during their decay (e.g. lOian et all Il999t iDomainko & Rufferu 
120051 120081) . The strength of these nuclea r lines depend o n 
the decay properties of the respective nuclei (lOian et all 1 1 999h . 
Promising nuclei in terms of detectability of the gamma-ray line 
emission comprise a half-life time r comparable to the age of 
the remnant. The estimated flux values in the table are given for 
a mass of the respective nuclei of 10~ 5 M , an age of the rem- 
nant of 10 4 , years and a distance to Terzan 5 of 5.9 kpc. Even 
for such an advantageous case, the line brightness would not ex- 
ceed 10~ 8 y ctrT 2 s _1 , and detection would be challenging for 
the next generation of hard X -ray / soft gamma-ray instruments 
(e.g. lKnodlseder et alll2009b . 



4. Subrelativistic ejecta 

4. 1 . Pressure driven remnant 

During the merger of compact binaries, a small fraction of the 
NS matter (10~ 4 - 0.1 M ) is dynamical ly ejected from the 
system with subrelativistic veloc i ties (e .g. IRuffert et all Il996t 
iLee & Kluzniakt Il999t iRosswogL 120051) . When assuming a ve- 
locity of the ejected material of c/3 (c speed of light), its ki- 
netic energy would be 10 49 - 10 52 ergs. After the ejecta has dis- 
placed a comparable mass of ambient medium (n as 0.1 cm -3 ), 
most of the k inetic energy is converted into thermal energy (e.g. 
iDorfil 1 19901) and the remnant expands farther in the pressure- 
driven phase dDomainko & Rufferll2005ll2008b . It should, how- 
ever, be noted that a large population of hadronic cosmic rays 
produced in the ultrarelativistic outflow may alter the evolu tion 
of the non-relativistic blast wave (e.g. IChugai et all 1201 lb . In 
the pressure-driven stage, the remnant is filled with hot ther- 
mal plasma that might be observable in thermal X-ray emission. 
Cooling of the thermal plasma appears to be unimportant in a 
remnant of the proposed age of a few times 10 4 years, since ma- 
terial with a density of 0.1 ctrT 1 and temperature of 10 7 (10 6 ) 
K would cool on a timescale of abo ut 10 7 (3 x 10 5 ) years when 
X-ray line emission is considered dSarazinl 1 19861) . Not detect- 
ing an additional thermal component above the galactic diffuse 
X-ray emi ssion and the exte nded nonthermal X-ray source with 
Chandra dEger et all l2010h would favor a content of hot ther- 
mal plasma in the potential remnant at the lower end of model 
predictions. A pressure driven remnant with energetics of 10 49 
ergs and an assumed age of 10 4 years would extend to a radius 
of about 10 pc corresponding to an angular radius of 6' at a 
distance of 5.9 kpc. It would be filled with thermal plasma of 
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E y 




nucleus 


[10 3 yr] 


[keV] 


[10~ 7 y cm -2 s- 1 ] 


22b Ra 


2.31 


609 


0.01 


229 Xh 


10.6 


40.0 


0.04 


243 Am 


7.37 


74.7 


0.09 



Table 1. Properties of gamma ray lines in a potential merger 
remnant in Terzan 5. 



5. Ionizing radiation 

5.1. Astrophysical signatures 

GRBs and their afterglows exhibit strong ionizing radiation 
which should alter the equilibrium state of the surrounding in- 
terstellar medium (ISM). Signatures of such an ionizing event 
should be obser vable for several 1 4 years in form of high- 
ionization lines dPerner et al 1 l2000b . Terzan 5 is, however, lo- 
cated in a direction close to the galactic center (1 3.8°, b 1.7°), 
hence affected by severe galactic absorption (Ay = 7.72, 
lOrtolani et allll996h . which may challenge the detection of such 
emission lines. 



5.2. Terrestrial signatures 

The ionizing radiation of a GRB that is beamed towards 
the Earth may a l so alt er the state of the Earth's atmosphere 
dFishman & Inanl [1988'), and Gal actic GRBs should produce a 
substantial amount of NO y there dMelott et all 12004b. This can 
be deposited as nitrates in the polar ice caps (e.g. Thomas et all 
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l2005tlEizak et all 120071) . which represe nt an archive of the at - 
mospheric conditions for 7.4 x 10 5 years (lAugustin et al.Ll2~004l) . 
Similar signatures of nearby, Gala c tic sup ernovae have already 
been claimed (e.g. iMotizuki et all l2009h . Here the possibility 
of such a signature originating from a GRB in Terzan 5 is ap- 
praised. Considerations in this section are less well founded than 
discussions in previous sections. An e valuation of different ter- 
restrial NO y signatures can be found in lHarfoot et all d2007l) . 

To assess NO y deposition from a GRB in Terzan 5, the so- 
lar proton event that occ urred during Augu s t 1972 and which 
is imprinted in polar ice dZeller et all 1 19861; IShea et all 120061) 
is used as a reference. This specific even t produced 3.6 x 10 3 3 
molecules NO y in the middle atmosphere dJackman et all 120051) . 
It is found that a fluence of about 10 5 ergs/cm 2 of ionizing radia- 
tion is needed for a comparable NO y en hancement if it is ado pted 
that 35 eV are necessary per ionization dGehrels et all 120031) and 
that each i onization releases 1.2 5 N atoms, which rapidly pro- 
duce NOy dJackman et al 1 12001 . A GRB at an adopted distance 
of Terzan 5 of 5.9 kpc would have to release an isotropic fluence 
comparable to typical short bursts of several 10 50 ergs in ioniz- 
ing radiation to produce enough nitrates to leave significant sig- 
natures in polar ice cores. Thus such an event could indeed have 
left traces in polar ice but only on the level of more frequent solar 
proton events. However, since GRBs are beamed with a beam- 
ing factor /b in the range of 1 «: f h l < 100, the probability that 
the prompt GRB emission has hit the Earth would only be a few 
percent. 

Along with generating nitrates by X-rays and soft gamma- 
rays, photons from the GRB with energies > 10 MeV would 
produce long-living radionucl ides like 10 Be ( half-life time 1.5 
Myr) in the Earth atmosphere dThorsettl.119951) . which will after 
rain-out also be deposited in polar ice caps. Such high-energy 
emission has indeed been observed from some short GRBs (e.g. 
lAbdo et all l2010bl) . In this framework it is interesting to note 
that two possibly cos mogenic 10 Be anomali es have been found 
in antarctic ice cores dRaisbeck et all Il987l) . However, the de- 
tectability of such signatures from a potential short GRB in 
Terzan 5 will again be limited by the beaming of the high-energy 
radiation. 

6. Summary 

The VHE gamma-ray source discovered in the vicinity of 
Terzan 5 appears to be compatible with being the remnant of 
a short GRB induced by the merger of two compact stars. An 
incontrovertible proof is missing, but also no contradictory facts 
have been found for this scenario. Further observations in the X- 
ray and gamma-ray regime may help to strengthen the evidence 
for a GRB remnant in Terzan 5. In particular X-ray observations 
could probe the presence of thermal plasma heated by subrela- 
tivistic merger ejecta. The identification of a GRB remnant in our 
own Galaxy would provide important constraints on the GRB 
rate in the very local universe and may also help to explore the 
effect of such a highly energetic event on the Earth. 
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